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I. INTRODUCTION AND OVERVIEW 



In the Standard Model (SM), the mechanism of CP 
violation in weak interactions arises from the joint ef- 
fect of three mixing angles and the single irreducible 
phase in the three-family Cabibbo-Kobayashi-Maskawa 
(CKM) quark- mixing matrix The unitarity of the 
CKM matrix V implies a set of relations among its el- 
ements, Vij, with i — u,c,t and j = d,s,b. In partic- 
ular, V ud V* b + V cd V* b + V td V* b = 0, which can be de- 
picted in the complex plane as a unitarity triangle whose 
sides and angles are related to the magnitudes and phases 
of the six elements of the first and third columns of 



*Now at the University of Tabuk, Tabuk 71491, Saudi Arabia 
TNow at European Organization for Nuclear Research (CERN), 
Geneva, Switzerland 

J Also with Universita di Perugia, Dipartimento di Fisica, Perugia, 
Italy 

SNow at the University of Huddcrsficld, Huddersfield HD1 3DH, 
UK 

^Deceased 

"Now at University of South Alabama, Mobile, Alabama 36688, 
USA 

ttAlso with Universita di Sassari, Sassari, Italy 



the matrix, Vid and Va,- The parameter 7, defined as 
arg[— V u dV* b /VcdVc b ], is one of the three angles of this 
triangle. From measurements of the sides and angles 
of the unitarity triangle from many decay processes, it 
is possible to overconstrain our knowledge of the CKM 
mechanism, probing dynamics beyond the SM In 
this context, the angle 7 is particularly relevant since it 
is the only CP-violating parameter that can be cleanly 
determined using tree- level B meson decays j^j. In spite 
of a decade of successful operation and experimental ef- 
forts by the B factory experiments, BABAR and Belle, 7 
is poorly known due to its large statistical uncertainty. 
Its precise determination is an important goal of present 
and future flavor physics experiments. 

Several methods have been pursued to extract 7 0- 
9]. Those using charged B meson decays into D^K^ 1 
and DK * ± final states, denoted generically as £)(*)i^(*) ± j 
yield low theoretical uncertainties since the decays in- 
volved do not receive contributions from penguin dia- 
grams (see Fig. [T]). This is a very important distinction 
from most other measurements of the angles. Here, the 
symbol indicates either a D° (D*°) or a D° {D*°) 
meson, and K** refers to if*(892) ± states. The meth- 
ods to measure 7 based on P ± — > £)(*) K^*^ decays 
rely on the interference between the CKM- and color- 
favored b — > cas and the suppressed b — > ucs ampli- 
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tudes, which arises when the D° from a B~ — > D°K~ 
decay [l(| (and similarly for the other related B decays) 
is reconstructed in a final state which can be produced 
also in the decay of a D° originating from B~ — > D a K~ 
(see FigQJ. The interference between the b — > cus and 
b — > ucs tree amplitudes results in observables that de- 
pend on their relative weak phase 7, on the magnitude 
ratio tb = \A(b — > ucs)/ A{b — > cus)\ and on the relative 
strong phase 5b between the two amplitudes. In the case 
of a nonzero weak phase 7 and a nonzero strong phase 
5b, the B~ and B + decay rates are different, a manifes- 
tation of direct CP violation. The hadronic parameters 
tb and 5b are not precisely known from theory, and may 
have different values for DK ± , D*K ± , and DK*^ final 
states. They can be measured directly from data by si- 
multaneously reconstructing several D decay final states. 



B 




B 




FIG. 1: Dominant Feynman diagrams for the decays B — > 
D°K~ (left) and B~ -> D°K~ (right). The left diagram 
proceeds via b —> cus transition, while the right diagram pro- 
ceeds via b —¥ ucs transition and is both CKM- and color- 
suppressed. 

The three main approaches employed by the B factory 
experiments are: 

• the Dalitz plot or Giri-Grossman-Soffer-Zupan 
(GGSZ) method, based on three-body, self- 
conjugate final states, such as K®tt + tt~ Q; 

• the Gronau-London-Wyler (GLW) method, based 
on decays to CP-eigenstate final states, such as 
K+K- and K°tt° Q; 

• the Atwood-Dunietz-Soni (ADS) method, based 
on D decays to doubly-Cabibbo-suppressed final 
states, such as D° -> K + ir~ [§]. 

To date, the GGSZ method has provided the highest 
statistical power in measuring 7. The other two meth- 
ods provide additional information that can further con- 
strain the hadronic parameters and thus allow for a more 
robust determination of 7. The primary issue with all 
these methods is the small product branching fraction 
of the decays involved, which range from 5 x 10~ 6 to 
5 x IIP 9 , and the small size of the interference, propor- 
tional to r B « c F \V cs V* b \/\V us Vl\ 0.1, where c F « 0.2 



is a color suppression factor [llMl3l |. Therefore a pre- 
cise determination of 7 requires a very large data sample 
and the combination of all available methods involving 
different D decay modes. 



Recently, Belle [lj] and LHCb [151 have presented 
the preliminary results of the combination of their mea- 
surements related to 7, yielding 7 to be (68^4) and 

(7lijg)°, respectively. Attempts to combine the results 
by BABAR, Belle, CDF, and LHCb have been performed 
by the CKMfitter and UTfit groups 0. Their most re- 
cent results are (66 ± 12)° and (72 ± 9)°, respectively 

The BABAR experiment [l6[ at the PEP-II asymmetric- 
energy e + e~ collider at SLAC has analyzed charged B 
decays into DK ± , D*K ± and DK*^ final states using 
the GGSZ [IM1, GLW [2^12 and ADS HUH meth- 
ods, providing a variety of measurements and constraints 
on 7. The results are based on a dataset collected at a 
center-of-mass energy equal to the mass of the Y(AS) res- 
onance, and about 10% of data collected 40 MeV below. 
We present herein the combination of published BABAR 
measurements using detailed information on correlations 
between parameters that we have not previously pub- 
lished. This combination represents the most complete 
study of the data sample collected by BaBar and bene- 
fits from the possibility to access and reanalyze the data 
sample (see Sec. |H] for detail). 

Other analyses related to 7 [H-H^l or 2/3 + 7 [H, EH 
have not been included, because the errors on the exper- 
imental measurements are too large. 



II. 



INPUT MEASUREMENTS 



In the GGSZ approach, where D mesons are re- 
constructed into the K®tt + tt~ and K®K + K~ final 
states dMl, the signal rates for B ± ->■ and 
B^ DK*^ decays are analyzed as a function of the 
position in the Dalitz plot of squared invariant masses 
m 2 _ = m 2 (K^h~), rn\_ — m 2 (K^h + ), where h is either 
a charged pion or kaon (h = tt,K). We assume no CP 
violation in the neutral D and K meson systems and ne- 
glect small D° — D° mixing effects [H, tH)]) leading to 
A(m 2 _,m 2 i _) = ^(m^rn^), where A (A) is the D° (D°) 
decay amplitude. In this case, the signal decay rates can 
be written as 34] 



T-i(*) /■ 2 2 \ 

1 ± \m_,m + ) oc 



1-4= 



|2 , J*) 



L~\A T \' 2 + 2\Re[z i * ) AiA T 



T s ± (m 2 _,m 2 + ) oc 



\A±\ 2 + r 2 s± \A T \ 2 + 2 Re[z s± 4i T ] , (1) 



with A± = A(m 2 ± , '1 



) and A ± is the complex conjugate 



of A±. The symbol A for B ± -> D*K ± accounts for the 
different CP parity of the D* when it is reconstructed 
into Dir° (A = +1) and D-f (A = —1) final states, as a 
consequence of the opposite CP eigenvalue of the ir° and 



the photon [35fl . Here, r B ± and r s ± are the magnitude 



ratios between the b — > ucs and b —> cus amplitudes for 
B ± -> D { -*^K ± and B ± -> DK* ± decays, respectively, 
and 5 B , 5 S are their relative strong phases. The analysis 
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extracts the CP-violating observables [1J 



,(*) - J*) 



x± 

X s ± 



■ (*) 



(2) 



denned as the suppressed-to-favored complex amplitude 



ratios z± 



s± 



Kr-ie*^*^), for 



B ± — > D ( -*' > K ± and B ± -> DK* ± decays, respectively. 
The hadronic parameter k is defined as 

iSa _ fA c (p)A u (p)e^dp 
Ke = — , 

^/jAl(p)dpjAl(p)d P 

where A c (p) and A u (p) are the magnitudes of the b — > 
and b — > ucs amplitudes as a function of the 



(3) 



— > 

DK^n^ 1 phase space position p, and 5(p) is their relative 
strong phase. This coherence factor, with < k < 1 in 
the most general case and k = 1 for two-body B decays, 
accounts for the interference between B^ — > DK*^ and 
other B ± — > DK^n^ 1 decays, as a consequence of the 
A'* ± natural width [ijj]- In our analysis, k has been 
fixed to 0.9 and a systematic uncertainty has been as- 
signed varying its value by ±0.1, as estimated using a 
Monte Carlo simulation based on Dalitz plot model of 
S ± — > DK^tt^ decays [18j. Thus, the parameter 8 S is an 
effective strong-phase difference averaged over the phase 
space. 



„(*) 



.,(*) 



TABLE I: CP-violating complex parameters z*j* = x±' +iy±' 
and z s ± = x 3 ± + iy a ±, measured using the GGSZ tech- 
nique [IJ]. The first uncertainty is statistical, the second is 
the experimental systematic uncertainty and the third is the 
systematic uncertainty associated with the D° decay ampli- 
tude models. The sample analyzed contains 468 million BB 
pairs. 





Real part (%) 


Imaginary part (%) 


z_ 


6.0 ±3.9 ±0.7 ±0.6 


6.2 ± 4.5 ±0.4 ±0.6 


z+ 


-10.3 ±3.7 ±0.6 ±0.7 


-2.1± 4.8 ±0.4 ±0.9 


z*_ 


-10.4 ±5.1 ±1.9 ±0.2 


-5.2± 6.3 ±0.9 ±0.7 


z + 


14.7 ± 5.3 ± 1.7 ±0.3 


-3.2± 7.7 ±0.8 ±0.6 


Z s - 


7.5 ±9.6 ±2.9 ±0.7 


12.7 ± 9.5 ±2.7 ±0.6 




-15.1 ±8.3 ±2.9 ±0.6 


4.5 ± 10.6 ±3.6 ±0.8 



12-dimensional (correlated) Gaussian probability density 

function (P.D.F.), defined in the space of the Zj^ and 
z s ± measurements from Tabled] After this combination 
has been performed, the values of 7 and of the hadronic 
parameters r p \ nr Sl 8g, and S s can be obtained. 

The D decay amplitudes A± have been determined 
from Dalitz plot analyses of tagged D° mesons from 
D* + — >• D°n + decays produced in e + e~ — )• cc events [H, 
|36| . assuming an empirical model to describe the varia- 
tion of the amplitude phase as a function of the Dalitz 
plot variables. A model independent, binned approach 
also exists 0, H3] , which optimally extracts information 
on 7 for higher statistics samples than the ones available. 
This type of analysis has been performed as a proof of 
principle by the Belle collaboration (38|, giving consis- 
tent results to the model-dependent approach 39]. The 
LHCb collaboration has also released results of a model- 
independent GGSZ analysis [40(. 

In order to determine 7 with the GLW method, the 
analyses measure the direct CP-violating partial decay 
rate asymmetries 



A 



(*) = 
cp± — 



T(B~ 



D { * } K- 



-r(s H 



nO) 

LJ CP± 1X 



) 



T(B~ D 



(*) 

CP± 



K~ 



D 



(*) 

CP± 



K+) 



(4) 



and the ratios of charge-averaged partial rates using D 
decays to CP and flavor eigenstates, 



R 



(*) 

CP± 



T(B- -> D%l ± K-)+T(B+ 



D 



(*) 

CP± 



K+) 



T(B- -> D(*)°K-) + T(B+ -> D(*)°K+) 



(5) 



Table U summarizes our experimental results for the 

(*) 

CP- violating parameters z ± and z s ±. Complete 12 x 12 
covariance matrices for statistical, experimental system- 
atic and amplitude model uncertainties are reported in 
Ref. [TtJ • The z± and z s ± observables are unbiased and 
have Gaussian behavior with small correlations, even for and 
low values of r B , nr s and relatively low statistics sam- 
ples. Furthermore, their uncertainties have minimal de- 
pendence on their central values and are free of physical 
bounds pjjj . These good statistical properties allow for 
easier combination of several measurements into a sin- 
gle result. For example, the rather complex experimen- 
tal GGSZ likelihood function can be parameterized by a 



where D^p ± refers to the CP eigenstates of the me- 
son system. We select D mesons in the CP-even eigen- 
states 7r~7r + and K~K + (Dcp+), in the CP-odd eigen- 
states K®it , K®(f), and K®uj (Dcp-), and in the non-CP 
eigenstate K-n+ (D° from B~ -)■ D°h.-) or K + n~ (D° 
from B + — » D°h + ). We recontruct D* mesons in the 
states Dtt° and Dj. The observables A S CP ± and R s cp± 
for B^ — > DK*^ decays are defined similarly. 

For later convenience, the GLW observables can be re- 
lated to z ± ' and z s ± (neglecting mixing and CP violation 
in neutral D decays) as 



.4 



(*) 

CP± 



±- 



(*) (*) 
x_ — Xt_ 



l + |z(*)| 2 ±(zL* j +<* J )' 



R { cL = l±|zW| 2 ±(x W ±. 



.(*) 



(G) 



(7) 



where |z'*'| 2 is the average value of |z+ | 2 and |z_ | 2 . For 
_B ± — > DK*^ decays, similar relations to Eqs. © and ([7]) 
hold, with k = 1, since the effects of the non-K* B — > 
DKtt events and the width of the K* are incorporated 
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into the systematic uncertainties of the A s CP:i _ and R S CP ± 
measurements [23 |. 

Table [TT| summarizes the results obtained for the GLW 
observables. In order to avoid overlaps with the sam- 
ples selected in the Dalitz plot analysis, the results for 
B 33 — > Dqp-K 33 decays are corrected removing the 
contribution from Dqp- — > K®(f>, <fi — > K + K~ candi- 
dates ED. For the decays B 33 -4 D* CP _[D C p-Tr°]K ± , 



B 3 



D. 



CP- 



[L> CP _7]i^ ± , and B 33 -> D CP -K 



*i 



such 



information is not available. In this case, the overlap 
is accounted for by increasing the uncertainties quoted 
in Refs. [2l|, [22| by 10% while keeping the central val- 
ues unchanged. The 10% increase in the experimen- 
tal uncertainties is approximately the change observed 
in B 33 — > Dqp-K 33 decays when excluding or including 
D — > K®<f) in the measurement. The impact on the com- 
bination has been found to be negligible. 



TABLE II: GLW observables measured for the B ± -> DK ± 



(based on_467 million BB pairs) 20], B ± ->• D* K ± (383 
million BB pairs) 0, and B ± DK*± (379 million BB 
pairs) [22l ] decays, corrected removing the contribution from 
Dcp- — ► K%<f>, <f> K + K~ candidates. The first uncertainty 
is statistical, the second is systematic. 



CP- even 



CP-odd 



Rcp± 
Acp± 
Rcp± 

A* 

A CP± 

Rcp± 

A 3 

A CP± 



1.18 ±0.09 ±0.05 
0.25 ±0.06 ±0.02 
1.31 ±0.13 ±0.04 
-0.11 ±0.09 ±0.01 
2.17 ±0.35 ±0.09 
0.09 ±0.13 ±0.06 



1.03 ± 0.09 ± 0.04 
-0.08 ±0.07 ±0.02 
1.10 ±0.13 ±0.04 
0.06 ±0.11 ±0.02 
1.03 ±0.30 ±0.14 
-0.23 ±0.23 ±0.08 



As in the case of the GGSZ observables, A^p , , A s cp , , 

R CP ±i and R^p± have Gaussian uncertainties near the 
best solution, with small statistical and systematic cor- 
relations, as given in Ref. [10] for B ± — > Dcp-K 33 de- 
cays. The GLW method has also been exploited by the 
Belle H3, CDF [H[, and LHCb collaborations @, with 
consistent results. 

In the ADS method, the D meson from the fa- 
vored -± ys amplitude is reconstructed in the doubly- 



Cabibbo-suppressed decay if" 1 



while the D from 



the —t lies suppressed amplitude is reconstructed in the 
favored decay K + tt~ [22, l23(. The product branch- 
ing fractions for these final states, which we denote 
as B~ — > [K+tt-] d K-, B~ -> [K + n-]* D K~, B~ -> 
[K + 7T~]p,K*~ , and their CP conjugates, are small (~ 
10 -7 ). However, the two interfering amplitudes are of 
the same order of magnitude, allowing for possible large 
CP asymmetries. We measure charge-specific ratios for 
B + and B~ decay rates to the ADS final states, which 
are defined as 



R 



(8) 



and similarly for R±, where the favored decays B 



-tt + } D K- 



B- 



[K- 



and B 



[K~ 

[K ir + ]p,K* serve as normalization so that many sys- 
tematic uncertainties cancel. The rates in Eq. (O depend 
on 7 and the B decay hadronic parameters. They are re- 
latcd to and z s ± through 



R$ = 



(*) ■ 

B± 



2Xr D 



(*) r (*) ■ r 

x± cos od — y± sin r> 



(9) 



where r D = \A{D° ->• K+n-)/A(D° -> K~>k + )\ and 
8d are the ratio between magnitudes of the suppressed 
and favored D decay amplitudes and their relative strong 
phase, respectively. As in Eq. ([1]), the symbol A for 
B 33 — > D*K ± decays accounts for the different CP parity 
of D* — > Dir° and D* — > D~f. The values of rp, and 8p 
are taken as external constraints in our analysis. As for 
the GLW method, the effects of other B ± -)• DK^n 33 
events, not going through K* 33 , and the K* 33 width, 
are incorporated in the systematic uncertainties on R±. 
Thus, similar relations hold for these observables with 
k = 1. 

The choice of the observables R± (and similarly for 
Pj_ and R±) rather than the original ADS observ- 
ables Pads = {R+ + and A AU s = (R- - 
P + )/2Pads @ is motivated by the fact that the set of 
variables (Pads ; Aads) is not well behaved since the un- 
certainty on Aads depends on the central value of Pads ; 
while P_|_ and R- are statistically independent observ- 
ables. Although systematic uncertainties are largely cor- 
related, the measurements of P+ and P_ are effectively 
uncorrelated since the total uncertainties are dominated 
by the statistical component. 

We have also reconstructed B 33 — > [K^n^n^pK 33 de- 
cays [24[ from which the observables P^ 71 ^ have been 
measured, which are related to the GGSZ observables as 



Rl- = r B± '< 



x [x± cos S K7r7r o - y± sin S Knw o] , (10) 

where n KlTn o is a D decay coherence factor similar to 
that defined in Eq. © for the B ± -> DK^n 33 decay, 
and where r^^-n and Sk-k-k are hadronic parameters for 
D° — > K^^n^ir decays analogous to rp, and 6p>. 

Table IIIII summarizes the measurements of the ADS 
charge-specific ratios for the different final states. Con- 
trary to the case of the GGSZ and GLW observables, 
R± , R±, and R± v do not have Gaussian behavior. 
The experimental likelihood function for each of the 
four decay modes, shown in Fig. [2] for P ± — > DK 32 
and B 30 — > D* K 33 decays, is well described around the 
best solution by an analytical P.D.F. composed of the 
sum of two asymmetric Gaussian functions. For the 
B 33 — > DK* 33 channel, we use instead a simple Gaus- 
sian approximation since in this case the experimental 
likelihood scans are not available. The effect of this ap- 
proximation has been verified to be negligible, given the 



small statistical weight of this sample in the combination. 
Measurements using the ADS technique have also been 
performed by the Belle S|, CDF 0, and LHCb 
collaborations [43| . with consistent results. 



TABLE III: ADS observables included into the combination 
for B ± ->■ DK ± with D ->■ A'tt (based on 467 million BB 
pairs) and D — > K-kty (based on 474 million BB pairs), 
B ± ->■ _D*J^ ± _(467 million BBpairs), and B± -» DK*^ 1 
(379 million BB pairs) decays [22 - [24l ]. The first uncertainty 
is statistical, the second is systematic. 



B i 



R± 0.022 ± 0.009 ± 0.003 

i?± [Dir°] 0.005 ± 0.008 ± 0.003 

R± [Dj] 0.009 ± 0.016 ± 0.007 

R s ± 0.076 ± 0.042 ± 0.0II 

D-fSTvrTr n nn r +0.012 +0.001 

n ± U.UUJ _ .010 -0.004 



0.002 ± 0.006 ± 0.002 
0.037 ± 0.018 ± 0.009 
0.019 ±0.023 ±0.012 
0.054 ±0.049 ±0.011 
n m 9 +0 012 +0.002 

u - ulz -0.010 -0.004 



III. 



OTHER MEASUREMENTS 



Similar analyses related to 7 measurement have 
been carried out using the decay B~ — > DK~ 



with the D — !> 
neutral B decay 



_o 



the 



B 



11 



K- 



r + 



with D 



hnal state [25|, and 
£>F*(892]° if* (892)° 

|| and D -i 
27]. For neutral B d& 



K°tt+it 



cays, tb is naively expected to be larger (« 0.3) be- 
cause both interfering amplitudes are color suppressed 
and thus cp ~ 1. However, the overall rate of events 
is smaller than for B~ -+ DK*~ decays. The flavor 
of the neutral B meson is tagged by the charge of the 
kaon produced in the K*° decay, K*(892)° -+ K~tt + or 
K*{892)° ->■ X+tt". 

Experimental analyses of the time-dependent decay 
rates of B -> D^T^i and £ -)• £>^£(770) ± decays have 
also been used to constrain 7 [28|, [29|. In these decays, 
the interference occurs between the favored b — ¥ cud and 
the suppressed b — > ucd tree amplitudes with and without 
B° — B° mixing, resulting in a total weak phase differ- 
ence 2/3 + 7 [20]; where (3 is the angle of the unitarity 
triangle defined as arg[~V c dV* b /V t dV t * b }. The magnitude 
ratios between the suppressed and favored amplitudes 
r D (,) v and rjy p are expected to be « 2%, and have to 
be estimated either by analyzing suppressed charged B 
decays (e.g., B + — > D + ir°) with an isospin assumption 
or from self-tagging neutral B decays to charmed-strange 
mesons (e.g., B° D+7r~) assuming SU(3) flavor sym- 
metry and neglecting contributions from IF-exchange di- 
agrams [30I ] . Performing a time-dependent Dalitz plot 
analysis of B — > D^K ^ decays [3l| could in principle 
avoid the problem of the smallness of r. In these decays 
the two interfering amplitudes are color suppressed, and 
tb is expected to be « 0.3 but the overall rate of events 
is too small with the current data sample. 




» 0.05 



FIG. 2: Experimental likelihoods as functions of the ADS 
charge-specific ratios R± (a,b), R± [B7r°] (c,d), R± [B7] (e,f), 
and R± 7Tn (g,h), from Refs. [23l . [24| . including systematic 
uncertainties. The P.D.F.s are normalized so that their max- 
imum values are equal to 1. These distributions are well pa- 
rameterized by sums of two asymmetric Gaussian functions 
with mean values as given in Table [TTTl 



In both cases, the errors on the experimental measure- 
ments are too large for a meaningful determination of 7, 
and have not been included in the combined determina- 
tion of 7 reported in this paper. However, these decay 
channels might provide important information in future 
experiments. 



IV. 



COMBINATION PROCEDURE 



We combine all the GGSZ, GLW, and ADS observables 
(34 in total) to extract 7 in two different stages. First, we 
extract the best-fit values for the CP-violating quantities 
whose definitions correspond to those for 



and z s ± 
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the quantities z± and z s ± of the GGSZ analysis given in 
Eq. ©. 

Their best-fit values are obtained by maximizing a 
combined likelihood function constructed as the product 
of partial likelihood P.D.F.s for GGSZ, GLW, and ADS 
measurements. The GGSZ likelihood function uses a 12- 
dimensional Gaussian P.D.F. with measurements z± and 
z s ± and their covariance matrices for statistical, exper- 
imental, and amplitude model uncertainties, and mean 
(expected) values and z s ±. Similarly, the GLW like- 
lihood is formed as the product of four-dimensional Gaus- 
sian P.D.F.s for each B decay with measurements A y cp± , 

A S CP ±, Rcp±> Rcp± an d their covariance matrices, and 
expected values given by Eqs. © and (J7]) after replacing 



the and z s ± observables by the z±' and z s ± parame- 
ters. Finally, the ADS P.D.F. is built from the product of 
experimental likelihoods shown in Fig. [5J With this con- 
struction, GGSZ, GLW, and ADS observables are taken 
as uncorrelated. Similarly, the individual measurements 
are considered uncorrelated as the experimental uncer- 
tainties are dominated by the statistical component. 

The combination requires external inputs for the 
D hadronic parameters td, <5d, r K7:7T o, S K7T7: o, and 
kktttt - We assume Gaussian P.D.F.s for rjj = 0.0575 ± 
0.0007 32] and r K7I7I « = 0.0469±0.0011 [13], while for the 
other three we adopt asymmetric Gaussian parameteri- 
zations based on the experimental likelihoods available 
either from world averages for 8 D = (202. Ot^) [32[ or 
from the CLEOc collaboration for S K7r7T o = (47±#)° and 
%„o = 0.84 ± 0.07 [11]. The values of S D and S K7r7V o 
have been corrected for a shift of 180° in the definition 
of the phases between Refs. [2J, [2J] and Refs. [32|, l48j . 
The correlations between rjj and Sp, and between KKmr 
and Skttk , are small and have been neglected. All five 
external observables are assumed to be uncorrelated with 
the rest of the input observables. 

The results for the combined CP- violating parameters 
z± and z s ± are summarized in Table HVl Figure [3] shows 
comparisons of two-dimensional regions corresponding to 
one-, two-, and three-standard-deviation regions in the 
z±, z?j_, and z s ± planes, including statistical and sys- 
tematic uncertainties, for GGSZ only, GGSZ and GLW 
methods combined, and the overall combination. These 
contours have been obtained using the likelihood ratio 
method, — 2Aln£ = s , where s is the number of stan- 
dard deviations, where 2 A In C represents the variation 
of the combined log-likelihood with respect to its maxi- 
mum value (47| . With this construction, the approximate 
confidence level (C.L.) in two dimensions for each pair 
of variables is 39.3%, 86.5%, and 98.9%. In these two- 
dimensional regions, the separation of the B~ and B + 
positions is equal to 2rs|sin7|, 2r^|sin7|, 2/«r s |sin7| 
and is a measurement of direct CP violation, while the 
angle between the lines connecting the B~ and B + cen- 
ters with the origin (0, 0) is equal to 27. Therefore, the 
net difference between x + and x_ observed in Table IIVI 
and Fig. [3J is clear evidence for direct CP violation in 



=(*) 



B ± DK± decays. 

In Fig. |31 we observe that when the information from 
the GLW measurements is included the constraints on 
the best fit values of the parameters are improved. How- 
ever, the constraints on y± are poor due to the quadratic 
dependence and the fact that < 1. This is the reason 
why the GLW method alone can hardly constrain 7. Sim- 
ilarly, Eq. © for the ADS method represents two circles 
in the (x±,y±) plane centered at {tb cos 5d,td sin Sd) 
and with radii v/ R±. It is not possible to determine 
7 with only ADS observables because the true (x±,y±) 
points are distributed over two circles [I9J]. Therefore, 
while the GLW and ADS methods alone can hardly de- 
termine 7, when combined with the GGSZ measurements 
they help to improve significantly the constraints on the 
CP- violating parameters z±, z* ± , and z s ±. 



TABLE IV: CP-violating complex parameters z± = x± + 

iy± and z s ± = x s ± + iy a± obtained from the combination 
of GGSZ, GLW, and ADS measurements. The first error is 
statistical (corresponding to — 2Aln£ = 1), the second is the 
experimental systematic uncertainty including the systematic 
uncertainty associated to the GGSZ decay amplitude models. 



Real part (%) 



Imaginary part (%) 



z_ 

z*L 
1+ 

Z s - 



8.1 ± 2.3 ± 0.7 
-9.3 ± 2.2 ± 0.3 
-7.0 ±3.6 ±1.1 
10.3 ± 2.9 ± 0.8 
13.3 ± 8.1 ± 2.6 
-9.8 ± 6.9 ± 1.2 



4.4 ± 3.4 ± 0.5 
-1.7 ±4.6 ±0.4 
-10.6 ±5.4 ±2.0 
-1.4 ±8.3 ±2.5 
13.9 ±8.8 ±3.6 
11.0 ± 11.0 ±6.1 



V. INTERPRETATION OF RESULTS 

In a second stage, we transform the combined (x± , y ± ) , 
(x±,y±), and (x s ±,y s± ) measurements into the physi- 
cally relevant quantities 7 and the set of hadronic pa- 
rameters u = (rB, 7-r, nr Sl 63, fi, S s ). We adopt a fre- 
quentist procedure 1501 to obtain one-dimensional con- 
fidence intervals of well-defined C.L. that takes into 
account non-Gaussian effects due to the nonlinearity 
of the relations between the observables and physical 
quantities. This procedure is identical to that used in 

Refs. [13, H mil nit 

We define a x 2 function as 



X 2 ( 7 ,u) = -2Aln£( 7 ,u) 

= -2[ln£( 7 ,u) - In £ n 



(11) 



where 2Aln£(7, u) is the variation of the combined 
log-likelihood with respect to its maximum value, with 
the z^ and z s ± expected values written in terms of 7 
and u, i.e., replacing z^ and z s ± by r^e 1 ^ >±7 ^ and 
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FIG. 3: (color online). Two-dimensional — 2Aln£ = s 2 contours (up to three standard deviations, i.e., s = 1,2,3) in the z± 
(left column), zj (center column), and z s ± (right column) planes, for the GGSZ measurement only (top row), the GGSZ and 
GLW combination (middle row), and the GGSZ, GLW, and ADS combination (bottom row). The solid (blue) and dashed (red) 
lines correspond to B~ and B + decays. 



nr s e l ^ Ss±1 \ respectively. To evaluate the C.L. of a cer- 
tain parameter (for example 7) at a given value (70), we 
consider the value of the x 2 function at the new mini- 
mum, x r 2 „i n (7o,uo), satisfying Ax 2 (70) = Xmm(7o,u ) - 
Xmin — 0- I R a purely Gaussian situation, the C.L. is 
given by the probability that Ax 2 (70) is exceeded for a 
X 2 distribution with one degree of freedom, 1 — C.L. = 
Prob[Ax 2 (7o); v = 1], where Prob[Ax 2 (7o); v = 1] is the 
corresponding cumulative distribution function (this ap- 
proach is later referred to as "Prob method") [47|. In 
a non-Gaussian situation one has to consider Ax 2 (70) 
as a test statistic, and rely on a Monte Carlo simula- 
tion to obtain its expected distribution. This Monte 
Carlo simulation is performed by generating more than 



10 9 samples (sets of the 39 GGSZ, GLW, ADS, and D 
decay observable values), using the combined likelihood 
evaluated at values (70, Uo), i.e., £(70, Uq). The confi- 
dence level C.L. is determined from the fraction of exper- 
iments for which Ax' 2 (70) > Ax 2 (70), where Ax' 2 (70) = 
x' 2 (7o,Uq) — Xmin f° r eacn simulated experiment is de- 
termined as in the case of the actual data sample. We 

adopt the Monte Carlo simulation method as baseline 

(*) 

to determine the C.L., and allow < r B ,nr s < 1 and 
-180° < 7,4* 3 A < 180 °- 

Figure 2] illustrates 1 — C.L. as a function of 7, rh\ 

nr s , Sg\ and S s , for each of the three B decay chan- 
nels separately and, in the case of 7, their combination. 
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The combination has the same twofold ambiguity in the 
weak and strong phases as that of the GGSZ method, 

( 7 ;4* ) ,(5 S ) -> (7 + 180°; + 180°, S s + 180°). From 
these distributions, we extract one- and two-standard- 
deviation intervals as the sets of values for which 1 — C.L. 
is greater than 31.73% and 4.55%, respectively, as sum- 
marized in Table |Vl When comparing these intervals to 
those obtained with the GGSZ method only, also shown 
in Table El we observe that the combination helps im- 
proving the constraints on r B and nr s , but not those on 
7. To assess the impact of the GLW and ADS observables 
in the determination of 7, we compare 1 — C.L. as a func- 
tion of r B and 7 for all B decay channels combined using 
the GGSZ method alone, the combination with the GLW 
measurements, and the global combination, as shown in 
Fig. [5] While the constraints on are clearly improved 
at the one- and two-standard-deviation level, and to a 
lesser extent on r B , their best (central) values move to- 
wards slightly lower values. Since the uncertainty on 7 
scales approximately as 1 jr B , the constraints on 7 at 
68.3% and 95.4% C.L. do not improve, in spite of the 
tighter constraints on the combined measurements shown 
in Fig. [3] However, adding GLW and ADS information 
reduces the confidence intervals for smaller 1 — C.L., as 
a consequence of the more Gaussian behavior when the 



significance of excluding r 



(*) 



increases. Thus, for 



example, in the region close to four standard deviations, 
the GGSZ method alone does not constrain 7, while the 
combination is able to exclude large regions. 



TABLE V: 68.3% and 95.5% 1-dimensional C.L. regions, 
equivalent to one- and two-standard-deviation intervals, for 
7, 5g , 5 S , rg, and nr s , including all sources of uncertainty, 
obtained from the combination of GGSZ, GLW, and ADS 
measurements. The combined results are compared to those 
obtained using the GGSZ measurements only, taken from 
Ref. [13] • The results for 7, S { * ] , and 5 3 are given modulo 
a 180° phase. 



Parameter 68.3% C.L. 95.5% C.L. 

Combination GGSZ Combination GGSZ 



7(°) 
r B (%) 

Tb (%) 

KT, (%) 

5 B n 
6 B n 
Ss n 



69] 



17 
16 

9.2±i:S 
1.9 
3.6 
+4.8 



10.6 

14.3 
105 



9.6 ±2.9 
,+4.2 
■3.9 

,+6.6 



4.!) 
hKi 

-66+1 
101 ± 43 



13.3 
14.9 

119_ 20 
-82 ± 21 
111 ±32 



(i.2 
-lit 



[41, 102] 
[6.0, 12.6] 
[3.0, 14.7] 
[3.3,25.1] 

[72, 139] 
-132, -26] 

[32, 166] 



[39, 98] 
[3.7, 15.5] 
[4.9,21.5] 
< 28.0 
[75, 157] 
-124,-38] 
[42, 178] 



The significance of direct CP violation is obtained by 
evaluating 1 — C.L. for the most probable CP conserving 
point, i.e., the set of hadronic parameters u with 7 = 0. 
Including statistical and systematic uncertainties, we ob- 
tain 1 - C.L. = 3.4 x 1CT 7 , 2.5 x 1CT 3 , and 3.6 x 1CT 2 , 
corresponding to 5.1, 3.0, and 2.1 standard deviations, 



o 




— B* ^ DK 1 

B* D*K* 

B* DK*± 
Combined 



150 -100 -50 




100 150 

Y (deg) 



o 




0.2 



0.3 



0.4 



0.5 



o 



1 

0.8 
0.6 
0.4 
0.2 



_ ' ' I ' ' ' ' I ' ' ' 


i i i i i i . i i i , i i i . i i . i i i . i _ 

^^^^ 











-150 -100 -50 



50 



100 150 



8^,5 s (deg) 



FIG. 4: (color online). 1 — C.L. distributions for the combi- 
nation of the GGSZ, GLW and ADS methods as a function of 
7 (top), Tg, and Kr s (middle), and Sg\ S B (bottom), includ- 
ing statistical and systematic uncertainties, for B ± — > DK ± , 
B ± — > D*K ± , and B ± — > DK* ± decays. The combination of 
all the B decay channels is also shown for 7. The dashed (dot- 
ted) horizontal line corresponds to the one- (two-) standard- 
deviation C.L.. 



for B ± M±, B ± D*K±, and B ± -)• DK** de- 
cays, respectively. For the combination of the three decay 
modes we obtain 1 — C.L. = 3.1 x 1CP 9 , corresponding 
to 5.9 standard deviations. For comparison, the corre- 
sponding significances with the GGSZ method alone are 
2.9, 2.8, 1.5, and 4.0 standard deviations [5lj], while with 
the GGSZ and GLW combination they are 4.8, 2.7, 1.8, 
and 5.4, respectively. 

The frequentist procedure used to obtain 7 and the 
hadronic parameters u is not guaranteed to have perfect 
coverage, especially for low values of r B \ r s . This is due 
to the treatment of nuisance parameters |50| . Instead 
of scanning the entire parameter space defined by 7 and 
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a = 0.679 ± 0.005 (0.955 ± 0.002) for the combination, 
and a = 0.670 ± 0.005 (0.950 ± 0.002) for the GGSZ 
method alone, for C.L. = 0.683 (0.954), respectively. 
For comparison purposes, the corresponding values using 
the Prob method are a = 0.641 ± 0.005 (0.941 ± 0.003) 
and a = 0.609 ± 0.005 (0.920 ± 0.003). While the Prob 
method tends to underestimate the confidence intervals, 
the Monte Carlo simulation method provides intervals 
with correct coverage, especially for the combination 
where the magnitude ratios between the suppressed and 
favored decays have more stringent constraints. 
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FIG. 5: (color online). Comparison of 1 — C.L. as a function of 
tb (top), r* B (middle), and 7 (bottom) for all B decay chan- 
nels combined with the GGSZ method only, the combination 
with the GLW measurements, and the global combination, 
including statistical and systematic uncertainties. Fhe hori- 
zontal lines represent the one-, two-, three- and four-standard- 
deviation C.L.. 



u (seven dimensions), we perform one-dimensional scans, 
in which, during MC generation, the nuisance parameters 
are set to their re-optimized best-fit values at each scan 
point. In order to evaluate the coverage properties of our 
procedure, we generate more than 10 9 samples with true 
values of (7, u) set to their best fit values, (7best, Ub os t), 
as given in Table [V] For each generated experiment, we 
determine 1 — C.L.' at 70 = 7best, as done previously with 
the actual data sample using the Monte Carlo simulation 
method. The statistical coverage a, defined as the prob- 
ability for the true value of 7 (70) to be inside the given 
1 — C.L. interval, is evaluated as the fraction of exper- 
iments with 1 — C.L.' larger than 1 — C.L.. We obtain 



In summary, using up to 474 x 10 6 BB decays recorded 
by the BABAR detector, we have presented a combined 
measurement of the CP-violating ratios between the 
-± ucs and -+ us amplitudes in processes B^ — > D^K^ 
and i? ± — > DK*^. The combination procedure max- 
imizes the information provided by the most sensi- 
tive 7 measurements and analysis techniques that ex- 
ploit a large number of D decay final states, includ- 
ing three-body self-conjugate, CP, and doubly- Cabibbo- 
suppressed states, resulting in the most precise measure- 
ment of these ratios. From the measurements of these 
ratios we determine 7 = (691^)° (modulo 180°), where 
the total uncertainty is dominated by the statistical com- 
ponent, with the experimental and amplitude model sys- 
tematic uncertainties amounting to ±4°. We also derive 
the most precise determinations of the magnitude ratios 
r B and nr s . The two-standard-deviation region for 7 is 
41° < 7 < 102°. The combined significance of 7 ^ 
is 1 — C.L. = 3.1 x 10~ 9 , corresponding to 5.9 standard 
deviations, meaning observation of direct CP violation in 
the measurement of 7. These results supersede our pre- 
vious constraints based on the GGSZ, GLW, and ADS 
analyses of charged B decays [i~7l - l2lT I23I [24j |. and are 
consistent with the range of values implied by other ex- 
periments [n-mms]. 
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